Impact of the transport supercurrent on the density of states in the weak link 
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The impact of the transport supercurrent on the quasiparticle density of states in the supercon- 
ducting thin film with the weak link is investigated. At the weak link the order parameter is locally 
suppressed due to the order parameter phase difference <j>. This results in the appearance of the 
zero-energy states, which are strongly influenced by the supercurrent especially at cj> close to it. The 
subgap density of states, dependent on both the supercurrent and the phase difference, is shown to 
modify the conductance characteristics of the structure. 

PACS numbers: 74.50.+r, 74.45.+C, 74.78.Na 

Properties of a mesoscopic superconductor are defined by its local density of states (LDoS). At the interface LDoS 
can be probed with the tunneling spectroscopy Particularly interesting are situations when there are states 
(so-called midgap states) within the bulk amplitude value of the order parameter Ao, that is at the quasiparticle 
energy s < Ao. This situation is ubiquitous in unconventional superconductors, where an interface results in the local 
suppression of the order parameter (due to its anisotropy) and in the appearance of the zero-energy states 2| ■ These 
zero-energy states are responsible for several interesting phenomena. First, when there is either external magnetic 
field or supercurrent injected, these midgap states create the countercurrent, which counterflows to the diamagnetic 
current in the former case, leading to the paramagnetic response 3, 4], or counterflows to the injected supercurrent 
in the latter case 5|. Second, when LDoS is probed with the tunneling spectroscopy, they result in the zero-bias 
conductance peak 3j, [6 1 . 

A weak link between two banks of conventional superconductors is similar to an interface of an unconventional 
superconductor in this context of local suppression of the order parameter 7] and of the creation of the midgap 



states 



8|. In recent years impact of the supercurrent on the density of states, and in particular, on the midgap 



states, has attracted the renewed attention (see e.g. Refs. S 



10 
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1J). In Refs. 



13[ we have studied first 



of the above mentioned effects in the weak link between current-carrying superconductors. And in this paper we 



"Electronic address: sshcvchcnko@ilt.kharkov.ua 



address the second problem of the appearance of the midgap states at the weak link between two current-carrying 
conventional superconductors. This subject is interesting both for further understanding of interference effects in 
mesoscopic superconducting structures and for possible applications for electronic devices, in which the current can 
be strongly controlled either electrodynamically (creating order parameter phase difference) or thermodynamically 
(injecting transport current) 



14| . In what follows we describe the model system to study the midgap states at the 
weak link between two current-carrying conventional superconductors. 

We consider the superconducting film with the impenetrable partition along y-axis, shown with the thick line in 
Fig. [TJ This partition (which can be e.g. a scratch in the film) has a break, which plays the role of the weak link in the 




FIG. 1: (Color online) Scheme for probing LDoS at the point contact between superconducting banks Sl and Sr with externally 
controlled order parameter phase difference <j) = 27r<I > e/ < I > o and externally injected supercurrent in parallel to the boundary 
between Sl and Sr. 



form of a slit between two superconducting banks Sl and Sr. The size of the weak link is assumed to be smaller than 
the coherence length (so-called pinhole model). The order parameter phase difference <j> at the weak link is created 
by the magnetic flux <I> e which pierces the loop connecting two banks. The magnetic flux can be created e.g. via the 
inductive coupling to a current-carrying coil. The supercurrent, tangential to the boundary between the banks Sl 
and Sr, is passed through the contacts, as in Ref. [12j. The LDoS is assumed to be probed by the scanning tunneling 
microscope's (STM) tip, placed over the weak link. Another contact for STM is not shown in Fig. [T]and we refer the 
reader to Ref. [12[ where both the current injection and the STM measurement as in Fig. [T]are discussed in detail. 

Thus, we consider a perfect point contact between two clean current-carrying singlet superconductors. The order 
parameter phase difference <fi is assumed to drop at the contact plane at x = 0. The homogeneous supercurrent 
flows in the banks of the contact along the y-axis, parallel to the boundary. The sample is assumed to be smaller 
than the London penetration depth so that the externally injected transport supercurrent can indeed be treated as 
homogeneous far from the weak link. Such a system can be quantitatively described by the Eilenberger equations. 



Taking transport supercurrent into account leads to the Doppler shift of the energy variable by Pfv s , where Pf is 
the Fermi momentum and v s is the superfluid velocity. The standard procedure of matching the solutions of the bulk 
Eilenberger equations at the boundary gives the Matsubara Green's function G™(u) at the contact at x — 5]. 
The analytic continuation (u> — > — ie + 7) of the component G^(w) = g( os, (uj) of G( )(w) gives the retarded Green's 
function, which defines LDoS at the boundary: 



N(e) = (N(s,p F )) = No (Re 5 (0) (e)) , (1) 

g^(e)=g ( -°\u J = -ie + 1 ), (2) 

(0) oj(tt L + fl R ) - isgn(v x )A L A R sm(j) 
g K [uj) — ^2 > \") 

flLfljl + OJ + AxAfl.COS</> 

here uj ~ TrT(2n + 1) are Matsubara frequencies, 7 is the relaxation rate in the excitation spectrum of the supercon- 
ductor, N(e, pf) is the angle-resolved DoS, iVo is the density of states at the Fermi level, (...) denotes averaging over 
the directions of Fermi momentum p^, Ax, ^ stands for the order parameter in the left (right) bank, uj = lj + ippv s , 



Ql.r. — yuj 2 + A| R . The poles of the retarded Green's function g^°\e) define the energy of the interface bound 
states. The direction-dependent Doppler shift pi?v s results in the modification of LDoS as it is discussed below. 

The method of the tunneling spectroscopy allows to access the LDoS by measuring the tunneling conductance G 
of the Superconductor - Insulating barrier - Normal metal (STM tip) structure. At zero temperature the dependence 
of the conductance on the bias voltage V is given by 1 1 

G{eV) = G N (D(p F )N(eV, p F )) , 

where Gn is the conductance in the normal state; D{p F ) is the angle-dependent Superconductor-Insulator-Normal 
metal barrier transmission probability. The barrier can be modelled e.g. as in Ref. [3| with the uniform probability 
within the acceptance cone |$| < i? c , where $ is the polar angle and the small value of i9 c describes the thick tunneling 
barrier. In this paper for simplicity we do not take into account the angle dependence of the transmission probability, 
assuming that the conductance is proportional to the LDoS: G(eV) tx N(e — eV). 

In the case of a weak link in the form of a constriction between two conventional superconductors with Ax = Ar = 
A = A (T,v s ) [13( Eq. © can be rewritten: 



a. 



where = Ao |cos(</>/2)| is the locally suppressed order parameter (so-called proximity gap) [7|. The retarded 
Green's function, determined with Eqs. (|2|) and |4}, gives for the energy of the Andreev bound states: ea = 
iA^ — pi?v s . Thus, there are the zero energy states, which are characterized by the values 4> &ud v s . The observation 
of the conductance characteristics of the system can be proposed as a test of the interface-induced transport-current- 



dependent quasiparticle states [3j, |8|, 



13] 



For the sake of generality we also write down below the Green's function for the current-carrying d-wave supercon- 
ductor. This can be easily done in the case of the specular reflection at the boundary, when the boundary between 
the current-carrying d-wave superconductor and the insulator can be modelled as the contact between two supercon- 
ductors with the order parameters given by A^ = A(#) and An = A(— <&) = A and with = 0. Then from Eq. (3) 
we have the following expression, which together with Eq. ([T]) describes the LDoS in the current-carrying d-wavc film 



as e.g. in Refs. 3 



n • n +s + a • A 



(5) 



where fl = + A 2 and Q 



u) 2 + A . This expression is valid for any relative angle \ between the a-axis 



and the normal to the boundary; in particular, at x = we have: </°'(a>) = uj/Q, and at \ — 7r /4 we have: 
g(°\u>) = n/2. Analogously Eq. (3) can be used to describe the LDoS at the boundary between two current-carrying 
d-wave superconductors as in Fig. 1 (see also in Ref. 5(). 

In what follows we consider the LDoS at the point contact between current-carrying conventional superconductors 
in details first analytically at 7 = and then numerically for 7^0. 

In the absence of impurities, that is in the limit 7 — > +0, we obtain the LDoS depending on the energy e and on 
the parameters Q = pfv s and <f>: 

N(e,Q,<j>) IV 



A 2 
^0 



Nn 



' (|e| - A ) + —As • S (\e\ - A^) ) , 



(6) 



\ 1 - A l 

where 8(...) and S(...) are the theta- and delta- functions; s = e — Pfv s ; A5 = A |sin(0/2)|. This after integration 
results in: 

N(s,Q,cf>) 



N 



^(±l) S5 n( £± )0(4-A 2 ) 



2Q 



(7) 



As ^ 

~2Q ^ 



(=Fl) sgn (e±) 6 (e 2 ± - A 2 ) arctan 



4-Ag n 



As +-9(Q-\e±AS 



Here e± = e ± Q: the first term describes the LDoS related to the continuum states, which coincides with the LDoS 
in the current-carrying homogeneous thin film 15 1, and the second term gives the contribution of the bound states. 



5 



Results of the numerical calculation at non-zero relaxation rate 7 are presented in Fig. 2. First, LDoS in the 




FIG. 2: (Color online) LDoS for the point contact between two current-carrying conventional superconductors for different 
values of the phase difference (j> and the dimensionless supercurrent velocity q = pfv 3 /Ao at zero temperature T = and 
non-zero relaxation rate 7/A0 =0.1. 



homogeneous current state (that is at (f> = 0) is shown in Fig. 2(a) to experience the suppression of the peak at 



15|. Then LDoS in 



e = Aq and gradual appearance of the midgap states with increasing the transport supercurrent 
the point contact without the supercurrent (q = 0) is plotted in Fig. 2(b). With increasing the phase difference <f> from 
to 7r the peak is shifted from e = A^ = o = Ao to e = A0 =7r = 0. Note that the peak at e = would appear as the 
so-called zero-bias conductance peak (ZBCP) in the STM measurements (see also in Ref. 16( about the appearance 
of ZBCP in point-contact Josephson junctions with ac biasing voltage). In Fig. 2(c) it is shown that the supercurrent 



121. 



results in the suppression and widening of this ZBCP similar to the ZBCP suppression in rf-wave superconductor 
Note that in concrete realizations the details of the modification of the ZBCP (widening, splitting) depend on several 
parameters, such as barrier anisotropy, surface roughness 



12 



171 ] . And finally in Fig. 2(d) we study the situation 



when the order parameter at the contact is significantly suppressed due to the phase difference at <f> close to n. In 
the absence of the supercurrent there are two midgap peaks (see also in Fig. 2(b)). With the gradual increase of 
the supercurrent these two peaks are shifted and modified so that first three-humped and then four-humped midgap 
structures appear, which are shown respectively for q = 0.4 and for q — 1. Similar humped structure of LDoS in 
homogeneous current-carrying superconductor (i.e. at <j> = 0) was studied in Ref. [llj. In the homogeneous situation, 



studied in Ref. 



ll|, the particularly interesting zero-energy states (which appear as ZBCP) take place when the 



6 

supercurrent is strong enough to depair electrons but smaller than the thermodynamic critical current, that is in 
rather narrow region: Ao < pfv s < 1.03Ao. We emphasize that in our case when the order parameter is locally 
suppressed (controlled) by the phase difference, the ZBCP appears at any value of the supercurrent at the phase 
difference defined by the relation = Ao cos(</>/2) < pfv s (see Figs. 2(b) and (d)). 

In conclusion, LDoS at the weak link between current-carrying superconductors have been studied. We propose this 
LDoS to be visualized with the STM. Two controlling mechanisms - external magnetic flux and injected supercurrent 
- allow to vary LDoS, and correspondingly the characteristics of the structure in different ways in the wide range of 
the parameters <p and v s . 
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